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done in biology. This gap is bridged by the work theorem. By using this theorem 
we show that the Barkhausen-type nonequilibrium noise in a process, repeated many 
times, can be combined to construct a special matrix S whose principal eigenvector 
provides the equilibrium distribution. For an interacting system S, and hence the 
equilibrium distribution, can be obtained from the free case without any requirement 
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1. Introduction 

A system in thermodynamic equilibrium has no memory of its past. Consequently 
there is no leading role for time in the ensemble based statistical mechanics except the 
subservient one to maintain equilibrium among the internal degrees of freedom and 
with external sources. This wisdom gets exploited in the dynamics based algorithms 
like Monte Carlo, molecular dynamics, stochastic quantization, to name a few, to attain 
equilibrium from any arbitrary state albeit in infinite time. Even a thermodynamic 
process involving changes in parameters is an infinite sequence of equilibrium states, 
and is therefore infinitely slow. A finite duration process, not destined to equilibrate 
at every instant of time, maintains a memory of the initial conditions or a short time 
correlation of states. The biased sampling of the phase space keeps these processes 
outside the realm of statistical mechanics and thermodynamics. In this equilibrium- 
nonequilibrium dichotomy, a work theorem [TJ |2l HI |5l [6] attempts to bridge the gap by 
providing a scheme for getting the thermodynamic free energy difference from a properly 
weighted nonequilibrium path integral jU [5] . 

We show in this paper that purely nonequilibrium measurements of work gives an 
operator S, defined on the phase or configuration space, whose normalized principal 
right eigenvector is the equilibrium probability distribution. Our result is valid for any 
number of parameters including temperature and interaction. With this extension we 
can get the equilibrium distribution by constructing a matrix S connecting any two 
allowed states of the system without any reference to equilibrium anywhere, thereby 
completely blurring the boundary between equilibrium and nonequilibrium. This finds 
direct application in out-of-equilibrium phenomena like hysteresis. 

Barkhausen noise is an example of nonequilibrium response of a ferromagnet as the 
magnetic field is changed at a given rate[8j M- By measuring the voltage induced in a 
secondary coil as the current in the primary coil wound around a ferromagnet is changed, 
one gets the time variation of the magnetization. The noisy signal one gets is not unique 
but stochastic in nature, reflecting the fluctuating microscopic response to the external 
field. Such signals have been analyzed in the past to extract information like avalanche 
statistics, material characteristics etc. Our results find a different use of the Barkhausen 
noise to construct the S matrix. Similar constructions for other cases like protein or 
DNA dynamics in vivo, pulling of polymers in single-molecule experiments, etc, call for 
new class of experiments to monitor the noise signals during these events. 

This paper is organized as follows: In Sec. El we recapitulate the work theorem, 
introduce the paths and discuss the connection between the work theorem and the 
histogram transformation of equilibrium statistical mechanics. In Sec. [3] we give a 
simple and general, dynamics independent proof of the relation between the equilibrium 
probability distribution and the work done in nonequilibrium paths. This relation in 
some form is already known [5j H] but our derivation allows us in generalizing the result 
to other cases involving temperature, interactions, etc. Sec. H] deals with the main result 
of this paper. There we prove the eigenvalue equation for S. A few examples are also 
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given there. How to get the operator S directly from experimental measurements of 
Barkhausen noise is also discussed here. Numerical verifications of some of the results 
are presented in Sec. [5] by taking the 2D Ising model as an example. We summarize in 
Sec. El 

2. Work theorem and path integral 

2.1. Work theorem 

Consider a classical system described by a Hamiltonian H(A,x) where A is an 
external field that couples to its conjugate, a microscopically defined quantity, x. The 
thermodynamic state is specified by temperature T and field A. Let us start with the 
system at A = in thermal equilibrium at temperature T. External field A is changed in 
some given way from to a final value A in a finite time r or in a finite number of steps 
n, letting the system evolve in contact with the heat reservoir. No attempt is made to 
ensure equilibrium during the process. The variation of x along the nonequilibrium path 
(x(t) vs t) and the instantaneous final (boundary of the path) value of x, Xb, when the 
field reaches A, are noted. The work done along a nonequilibrium path by the external 
source (as in ref. [2]) is 



in time r, and it varies from path to path. The difference between two definitions of 
work in the context of work theorem, one used in ref. |T] and the other in ref. [2], is 
discussed in ref. [3]. For the sake of notational simplicity we choose, 



where H is the energy for A = 0. There is not much loss of generality in choosing 
the form of Eq. [2] because A and x refer to any pair of conjugate variables so that 
x itself need not be a linear function of the internal coordinates. As an example, in 
an interacting spin problem in a magnetic field h (= A), H = Hq — h^2 k Sk where Sk 
is the spin variable at a site denoted by k, with x = J^ fc s&. Often A can be taken 
as the switching parameter to turn on a perturbation or interaction in a Hamiltonian 
H = H -H' with H A = H - AH'. 

The work theorem [H [2] provides the equilibrium free energy difference AF between 
the two states with A = and A = A, both at inverse temperature /3 = l/ksT (ks is 
the Boltzmann constant), from the nonequilibrium work done as 




(1) 



H = H + H^A, x) = H -Ax, 



(2) 



AF = ~]n(e- pw ), 



(3) 



where (...) denotes the average over all possible paths. 
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2.2. Paths: equilibrium and nonequilibrium 

We are using here a description of a state by the intensive parameters which actually 
characterize the surroundings. In equilibrium any system is expected to have the values 
of the intensive parameters same as that of the environment. A change in any of the 
parameters, say A, from Ao to A, would require heat and/or energy transfer. The work 
done on or by the system is determined by the change in the free energies, independent 
of the path of variation of the intensive parameters. This is expressed as 



where AF = F(/3, A) — F(f3, Ao). Here x eq (A) = J x P\(x)dx is the equilibrium average 
at the instantaneous values of the intensive parameters and P\(x) is the corresponding 
equilibrium probability distribution of x. This follows from the identification of the 
equilibrium value of x as x oq = —dF/ dA, in contrast to the conjugate ensemble definition 
A = dJ- ' jdx where T{J5, x) is the fixed-x ensemble free energy. 

For convenience, let us discretize the integrals. For example, for A G [Ao,A], we 
have a sequence (A , A 1; ...A n = A) and the continuum is recovered by taking the usual 
limit of n — > oo with max{AAj = A»+i — Aj} — > 0. The work done can be rewritten as 



By interchanging the sums over x and A, we define (i) a sequence {xi\i = 0, ...n} as 
instantaneous values, and (ii) a sequence- dependent work done as W = Ylii x i AAj, to 
reinterpret Eq. <^ as an average over these Xj's. Therefore, 



where V{xi} = JJ i P^(xi) is the joint probability of getting the particular {xi} sequence, 
because, for a thermodynamic process, there is no memory. Going over to the continuum 
limit, the thermodynamic process of varying A is now seen as equivalent to choosing 
a path in the configuration space and re-weight the paths according to the probability 
of its occurrence in the A-ensemble. The relation between the free energy change and 
work, Eq. (j4]), now gets a path integral meaning where the process takes the system 
over the microstates and one averages the work over individual paths. 

This thermodynamic connection is valid only in equilibrium. The work theorem 
generalizes this idea by replacing V{xi} by the nonequilibrium probability of getting a 
path and asserting 




(4) 




(5) 




(6) 



{Xi} 




(7) 



where J T>X stands for the normalized sum over paths, i.e., sum over intermediate x's 
with appropriate probabilities. 
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2. 3. Histogram transformation and infinitely fast process 

There is a fundamental transformation rule obeyed by the partition function, often used 
in numerical simulations as the histogram method[7]. This transformation connects the 
equilibrium probability distributions at two parameter values, A = Ao and A = A as 

P x (x) e P( x - x o)% 

where the sum in the denominator is over the allowed values of x. The denominator 
of the right hand side of EqJUis Z\/Z\ where Z\ is the partition function at inverse 
temperature (3, 

Z A = e-^e?**. (9) 

states 

From Eq. [TJ (A — \q)x can be taken as the work done in an instantaneous process that 
changes A from A to A without changing x. The probability of getting x for equilibrium 
at Ao is P\ (x) and therefore the sum in the denominator of Eq. [S] is the path integral 
of Eq. [7J because x does not change. This gives the work theorem. 

3. Equilibrium probability distribution 

We in this section use the discrete version of the process to re-derive the equilibrium 
probability distribution from the work theorem in a general and dynamics independent 
way. For the kind of nonequilibrium processes mentioned in Sec. 12.21 the equilibrium 
probability distribution of x at a parameter value A can be obtained from a weighted 
path integral [H [5] 

x{ 1 JVXe-» w ' 1 ' 

where x^ is the instantaneous boundary value at the end of the path, and the 
denominator is same as r.h.s. of Eq. [7J This is in the form of a path integral where the 
paths are weighted by a Boltzmann-like factor exp(— (3W). The same was established 
previously in specific cases like, the Master equation approach[2], the Feynman-Kac 
formula [5] and Monte Carlo dynamics [I]. 

The equilibrium average x eq is defined as 

1B .^ = ^(p^YU^-^), (id 



OA s^o \ Z J 5 \Z Z j 

where work theorem is to be used for the partition functions. 

The system starts in equilibrium at temperature T and A = 0, and then A is built 
up at constant T as a sequence of infinitely fast jump of AA = X/n, each jump followed 
by a finite time evolution in contact with the heat bath. Consider now two n-step 
processes, one process with final field A and another one with A — 5 (5 — > at the end). 
In fact, the second process is just a copy (replica) of the first one in every respect except 
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i=n—l A=X—AX 



i=l , A=AA. 
i=0, A=0 



replica 1 




replica 2 



Figure 1. Schematic representation of two replicas of same paths, each starting from 
A = and ending at A = A in replica 1 and at A = A — S in replica 2. Label i denotes 
the step number as A is changed in steps of X/n. Lines of different styles (dashed, 
dotted etc) represent different realizations of paths starting from different values of x. 
The vertical portion of a path is an instantaneous process (no change in x) and the 
horizontal part is under interaction with the surrounding (x evolves at a constant A) . 
Identically shaded lines in the two replicas have the same evolution. 



at the last stage (Fig. [I]). For the last jump, the change in A for replica 1 is A A while 
for replica 2 it is A A — 5. 

A path is specified or defined by the sequence {x t \ i = 0, ...n — 1}. The changes 
in Xi at any step is because of internal dynamics or exchange of heat with the external 
reservoirs. We do not need to let the system evolve once the field reaches the final 
desired value. Therefore, the sequence {xi \ i = 0...n — 1} is the same for both the 
replicas. The work done Wi,W2 along an n-step nonequilibrium path for replicas 1,2 
are related via 

W 2 = W 1 + 6x n - 1 , (12) 

with W\ is of the form given above Eq. ([6]). The work theorem of Eq. ([7]) when used in 
Eq. [TT] yields 

x cq = hm 



BY u e^fco AA < Xi 5 

r L^i paths ° 

jVX x h e~ pw 

= $VX e~P w ' ( Xb = Xn -^- ( 13 ) 

This shows that the equilibrium average can be expressed in terms of the boundary 
value with proper weightage of the paths. The above proof can be generalized to any 
moments of x. 

Now if V{x) is the distribution of x^, that gives the average in Eq. [13] 



x eq = (x) = J xV(x) dx, (14) 
then V(x) can be written as 
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as quoted in Eq. [TU] We now invoke the moment theorem[T2] which, in our case, 
states that for a probability distribution without sufficiently long tails, the moments 
uniquely specify the distribution. Since these conditions are satisfied by the equilibrium 
probability distributions for any finite system, the moment theorem applies. Since the 
moments from the nonequilibrium path integral are the equilibrium moments, V(x) is 
the equilibrium distribution: V(x) = P\(x). This completes the proof. 

3.1. Generalization 

In general, for a Hamiltonian of the form H = H({A a }, {X a }), the equilibrium 
distribution, P(E, x\, x 2 , ■■■), at some given parameter values, {A Q } and temperature 
can be obtained in the same way provided the paths start from an equilibrium 
state for H = H , where H gives the energy for all A a = and W is the total work 
done on the system along a nonequilibrium path, by each of the externally controlled 
parameters. E here corresponds to the energy from Hq only. Our starting Hq may 
be a free Hamiltonian for a mechanical system and can as well be zero for interacting 
spin-like systems. 

Consider the Hamiltonian H = ^Hq for a spin-like system (i.e. without any kinetic 
energy). In this case one of the {A a } could be the strength of interaction. Let's start 
with 7 = 0, i.e. the starting point is any random configuration of the free system or a 
non-interacting system, and then change 7 in some given way from 7 = to 7 = 1. We 
thus generate the equilibrium distribution of Hq at a particular /3, by doing a similar 
nonequilibrium path averaging. Note that everywhere we need the product j3W. So, 
we can discretize temperature instead of A and the process can be reinterpreted as 
cooling down to a finite temperature from an initial infinite temperature. In the usual 
formulation of work theorem, A refers to mechanical parameters such as the pulling 
force in AFM, which are under direct control of the experimentalists. In contrast, 
other intensive parameters such as temperature may not be controlled with this level of 
precision in experiments. But this finds various applications in numerical experiments. 
Such thermal quenches are quite common in numerical simulations and our results show 
how these can be harnessed to extract equilibrium information as well. The ensemble of 
states obtained in the above discussed way at the end of the path is not a representative 
sample of the equilibrium ensemble at the concerned temperature and field. However, 
the history-averaged distribution is the equilibrium distribution. The boundary states 
would relax to reach equilibrium via energy transfer to the reservoirs but that part of 
the process is not required. This difference becomes important and visible in systems 
exhibiting hysteresis as e.g. for a ferromagnet. 

3.2. Application to ferromagnet to get equilibrium magnetization curve 

The above-mentioned scheme can be used to get the equilibrium probability distribution 
or thermodynamic quantity from a process which is arbitrarily away from equilibrium 
and at all temperatures including phase transition points. Now we apply our result to 
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the case of hysteresis of a ferromagnet below the critical temperature (Tq). Consider a 
Hamiltonian: H = H — hM. The external magnetic field is varied from — ho to +h in 
a fixed manner and then reversed. (M) is calculated using Eq. [13] Below the critical 
temperature, magnetization (M) vs. magnetic field (h) curve shows a discontinuity at 
h = for infinite system size. For a finite system there is no discontinuity, M-h curve is 
continuous passing through the origin, and the slope of M-h curve at h — increases as 
system size increases. But, in reality, when experiments or simulations are done, instead 
of single retraceable curve passing through the origin we get a loop called hysteresis loop, 
no matter how slowly we vary the magnetic field. The common technique known to get 
the equilibrium curve is to connect the vertices of the sub-loops [9]. Here the weighted 
nonequilibrium path integral scheme is a way out to get the equilibrium magnetization 
curve. We verify this for Ising ferromagnet and discuss the observations about it in Sec. 

El 



4. Equilibrium probability distribution from an eigenvalue equation: 
Operator S 

In this section we derive the main result of this paper: equilibrium probability 
distribution as an eigenfunction of a nonequilibrium operator S. 
Using the discrete notation, we can write Eq. [10] as 

Px{x) = f^£ e-? w 5 Xb , x , (16) 

paths 

by using the work theorem, Eq. [3] that 

E ^ w = f^- ( 17 ) 

paths ^° 

Again, writing ^ paths = Yl Xi P*o( x i) X) pa ths> wnere the primed summation denotes the 
sum for fixed initial value of x = x- x with appropriate probability and P\ (x\) denotes 
the equilibrium distribution of X\ for A = A , we get, 

Px{x) = 7 fJ^Y. ' P ^0 e-? w S Xh , x . (18) 

x; paths 

Use the transformation rule for the partition function (Sec. I2.3p . 

!^ = £ P ^)e^)*, (19) 

to absorb Z\ Q /Z\ into the probability distribution. This transforms P\ (xi) into P\(xi), 
in Eq. [18] as 

= E E ' e"^-^-^ 6 Xh , x P x ( Xi ) (20) 

Xi paths 

= e^^-)- ( 21 ) 



■>:■) 



5P a = Pa, (22) 
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with F\ as a column vector of {P\(x)} and the matrix elements of S as 



C = ' e -/3W-^(A-A )xi 



(23) 



paths 



The summation in Eq. [23] is over all paths that start from an equilibrium distribution 
of A = Ao with value of x as x- x and end in a state with A = A and x — Xf , with proper 
normalization (denoted by prime). 

Although we use the simple Hamiltonian: H = H — A x in the construction, Eqf23l 
can be generalized for a Hamiltonian H = H + Hi(A, x), because Eq. [19] has the general 
form, 



Now we address the remaining problem - the normalization of the primed 
summation over paths in Eq. [23] This problem is inherited from Eq. [TT1 Note that 
the l.h.s. of Eq. [T7] should add up to 1 for A = A with W = 0. So we choose the 
hidden factor a posteriori by demanding proper normalization of the final probability 
distribution. This condition can be ensured in a process- or system-independent way by 
choosing $x,xi — f( x i) — 1> (Eql2"TT). i.e. by making the column sum of S independent 
of x\. By this normalization of the sum of each column to unity it is also guaranteed that 
the principal eigenvalue is 1. The corresponding right principal eigenvector has all the 
elements real and non-negative - a necessary condition to be a probability distribution 
and when normalized, such that sum of all elements is unity, this eigenvector gives the 
equilibrium probability distribution. 

The number of rows and columns in S is determined by the number of allowed 
values of x. For continuum of states, the matrix equation is to be replaced by an 
integral eigenvalue equation. 

Hence, in brief, the scheme to get the equilibrium distribution at some parameter 
value A and temperature (3~ l is as follows: Pre-fix some arbitrary or convenient-to-start- 
with initial parameter value Ao which will be same for all paths/experiments. Choose 
a microstate from the equilibrium distribution at field Ao and call its value of i as i;. 
Change the parameter value from Ao to A in some predetermined way and measure 
the work done by the external parameter on the system according to Eq. [1] Repeat 
the experiments several times and construct the matrix S using Eq. [23] Next, each 
column of the matrix is normalized to unity. The normalized principal eigen- vector is 
the equilibrium probability distribution, P\(x), at the field A. 

Eq. [22] is the main result of this paper and it is not restricted to one external 
parameter only and can be generalized to any parameter as mentioned above. The 
matrix S connects any two allowed states of the system without any reference to 
equilibrium anywhere and yet its principal eigen-vector determines the equilibrium 
distribution. Despite resemblance, there is no similarity either with the stochastic matrix 
of a Markov process or the adiabatic switching on of interaction in a quantum system 
because S is constructed out of a finite process and needs global information about the 
work done. 
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Another issue that comes up in this approach via S, is the question of ergodicity 
which connects the Gibbsian statistical mechanics with equilibrium thermodynamics. 
The nonequilibrium dynamics used to construct S may not respect ergodicity but the 
starting points for the paths in principle span the whole phase space, even in the case 
when one starts with a free non-interacting system. It seems ergodicity of the free 
noninteracting system is sufficient to generate the equilibrium distribution. 

4-1. Examples 

4-1.1. Example 1: Extreme cases Consider an extreme case: a completely equilibrium 
evolution of the system, where at each step the system reaches its equilibrium. Take 
a simple system: a single spin problem in magnetic field h and temperature /3 _1 : 
j3H = —Ks, where s = ±1 and K = (3h. For an n-step process, K varies from 
to nk in steps of k, and the column normalized S matrix can be calculated exactly 
where at each step the spin reaches the corresponding equilibrium state, as 



where P n &(±) is the equilibrium probability of finding ±1 spin at the n-th step. Thus 
for a completely equilibrium evolution of the system the elements of the matrix S are 
unique and, therefore, S has only one and unique eigenvector. In that case principal 
eigenvalue is 1 and all other eigenvalues are zero. We may conclude that a complete 
reducibility of S is the signature of a thermodynamic process. 

Eq. [24] is to be compared with the extreme nonequilibrium process as embodied 
in Eq. [HJ For this instantaneous change in A, S = I, the identity matrix, with no zero 
eigenvalues. 

If at each of these n steps, the system evolves for a time At in contact with the 
bath, then <S nj At S eq as At — > 00. The smallness of the rest of the eigenvalues would 
indicate how close to equilibrium the system is. 

The dynamics of a many body system might be compartmentalized into slow modes 
and fast modes, where the fast modes would equilibrate much more quickly than slow 
ones. How many such fast modes have actually equilibrated, can be gauged by the 
number of zero eigenvalues. The S matrix is not necessarily symmetric, though real and 
there is a possibility of pairs of complex conjugate eigenvalues, with their magnitudes 
going to zero as equilibrium is reached. 

4-1.2. Example 2: Barkhausen noise and matrix S We now show the practical 
feasibility of the operator method for a magnet by using the Barkhausen noise [H |9] 
as recorded through the output voltage across a secondary coil wound around a 
ferromagnetic material. Though Barkhausen noise has seen many applications, its use 
for equilibrium properties has not been anticipated. 
Consider the Hamiltonian 




(24) 



hM. 



(25) 
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Here magnetic field h and magnetization M correspond to A and x respectively. The 
field is varied from h\ to hf in a time interval r at a constant rate h. The Barkhausen 
effect is a noisy signal proportional to the change in magnetization, rj(t) = dM d f^ ■ So by 
integrating the Barkhausen noise up to time t one gets the nonequilibrium instantaneous 
magnetization of the material. Therefore, we can write the work related exponent in 
Eq. |23| as 



W + [h(r) - h(0)} M 



dt [ T]{t') dt', 
'o Jo 



(26) 



which, in a discretized form, looks like 
W+[hf- h] Mi 



n-l j 

-Afc£I>> (27) 

3=1 k=l 

where the Barkhausen noise at k-th step is r] k = M k — M^-\. Hence the matrix elements 
ShdiMi takes the form 

n-l j 

'exp \/3AhJ2J2^\ ' ( 28 ) 

expts. L 3=1 k=l 

expressed entirely in terms of the Barkhausen noise along the nonequilibrium paths. 
The primed summation over paths that start with M; and end at Mf includes proper 
normalization as mentioned earlier. 

To go to other cases, e.g., for the case of a polymer pulled at a constant rate of 
change of force, one needs to monitor the time variation of the pulled point displacement 
dx/dt vs t. This information can then be used in Eq. [28] to get the corresponding S. 



5. Numerical verification of results 



Our claims about the probability have been verified for the case of 2D Ising model on a 
square lattice, Lx L, where L is the size of the lattice with periodic boundary condition. 
Consider the Hamiltonian 

H = -jJ2 s kSl -hJ2s k , (29) 

<k,l> k 

where J is the interaction strength, h is the external magnetic field and s k = ±1 is the 
spin at fc-th site of a square lattice. Here ^2 <k ;> denotes the sum over nearest neighbor 
spins. Here J and h play the roles of external parameter (A) and ^2 <k t> s k si and J2k s k 
are the internal variables (x). 

We find equilibrium probability distribution for given J and h using weighted 
nonequilibrium path integral, normalizing the eigenfunction of S and compare those with 
the equilibrium probability distribution obtained from a usual Monte Carlo procedure. 
The overlap of the two distributions is determined by the Bhattacharyya coefficient |10j 
defined as 

BC = y/Ph(E,M)P eq (E,M) = 1 - e, (30) 

E,M 

with BC = for no overlap and BC = 1 for complete overlap. 



Thermodynamics as a nonequilibrium path integral 



12 



5.1. Numerical verification of the equilibrium probability distribution starting from a 
uniform distribution 
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Figure 2. Plot of the weighted distribution (a) Pj,h{E) vs. E and (b) Pj,h(M) vs. M 
(dotted line with circles) for varying J and ft with n = 20 and equilibrium distributions 
P eq (E) and P eq (M) (crosses) with J=l, ft. = 1 and (3 = 0.2 for a 8 x 8 lattice, showing 
that Pj th (E) = P eq {E) and Pj, h {M) = P eq (M). 

Let us take an 8 x 8 lattice and start from H = 0. Each time we start from a state 
chosen from a uniform distribution and reach the final state with J = 1 and /i = 1 in 
n-steps. At each i-th step, J is switched from J, to Jj + i and the external magnetic field 
from hi to h i+ i, 

A J = Jj + i — Ji — J/n and A/ij = — hi = h/n; 

keeping the spin configuration unchanged, and the amount of work done on the system 

Wi = -AJtEi - Ah t M t , 

is calculated where Mj is the magnetization and E { is s k$i at the i-th step. Then we 
let the system relax at that field hi, Ji and j3 for a while, but do not equilibrate. Thus 
the work along a path consisting of n steps is 

n-l 

W = -^2AJ i E i + Ah i M i , 

which is different for different paths. We find the weighted distribution 

/ VX e-P w 5(E h - E)6(M h - M) 



Pj,h(E,M) 



-pw 



(31) 



and then 



and 



Pj, h (M) = Y J P(E,M) 

E 

P Jth (E) = J2P(E,M). 

M 

It is observed that these distributions merge well with the corresponding equilibrium 
distributions and for Pj, h (E) (FigJ2(a)) and Pj, h (M) (Figfj^b)) we get e ~ 10~ 3 (Eq. 

Em. 
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5.2. Equilibrium magnetization curve using nonequilibrium path integral 

For this case lattice size is 8 x 8 and the interaction strength is kept fixed at J = 1. 
Each time we start from an equilibrium distribution of h = —ho. The field is varied 
from —ho to +ho in n steps. W(n) vs. n data are recorded and (M)(h) is calculated 
using Eq. [131 

We plot the weight averaged magnetization curve, (M)(h), along with the hysteresis 
loop, average magnetization over samples, against h for ho = 0.2 in Figj3]and ho = 2 in 

Fig. m 

60 
40 
20 

^ Q 

-20 
-40 
-60 

-0.2 -0.1 0.1 

h 

Figure 3. Plot of weighted average M(h) vs. h (black solid line) and hysteresis loop, 
simple averaged M vs. h (green dashed and blue dash-dotted lines) for a 8 x 8 lattice. 
The magnetic field h varies from —0.2 to +0.2 in 100 steps. Inset shows the hysteresis 
loop for the small (green and blues lines) field with respect to the large field (red double 
dash-dotted line) and the weight averaged magnetization for small field. 

A retraceable equilibrium curve is obtained as expected though the nominally 
averaged magnetization neither changes sign nor makes a complete loop (FigJ3l) [TT] . 
This reflects the fact that though in majority the magnetization does not reach the 
correct value, there are a few rare samples for which the spins do flip and these rare 
configurations, which are close to equilibrium, get more weight in the weighted path 
integral to give the correct equilibrium curve. 

For the larger field, we obtain a curve which is much narrower than the hysteresis 
curve (Fig|4j). The equilibrium curve obtained this way is still not a single curve. The 
width of the loop might be connected to the droplet time scale, and signals the need for 
a more careful sum over paths to take care of droplet fluctuations. 

5.3. Numerical verification of the eigenvalue equation 

We start from an equilibrium ensemble at inverse temperature (3 = 0.2 (kept fixed 
throughout the experiment), J = 1 and h = 0. Each time we start from a state chosen 
from its equilibrium distribution and reach the final state with J = 1 and h — 1 in n- 
steps in the same way described above and calculate the amount of work on the system 
at z'-th step: W{ = —AhiMi. We find the matrix elements: 

S MlM = £ ' e-W-PW* 5 MhiMf . (32) 

paths 
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Figure 4. Plot of weighted average M(h) vs. h (dashed lines) and hysteresis loop, 
simple averaged M vs. The magnetic field h (red solid line) for a 8 x 8 lattice, h varies 
from —2 to +2 in 100 steps. 



After the matrix is constructed, we normalize sum of each column to unity and 
find the normalized principal eigen-vector corresponding to the Principal eigenvalue 
1, which is guaranteed. We compare the normalized eigenfunction with the actual 
equilibrium distribution for L = 4 and 8. We see that these distributions merge with 
the corresponding equilibrium distributions for L = 4 (FigJS^a)) and L = 8 (Figj5](b)) 
with e ~ 1(T 4 (Eq. [30]). 
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Figure 5. Plot of the equilibrium distribution P eq (M) vs. M (boxes with dotted line) 
and normalized principal eigen-vector Ph(M) (dashed line with circles) with J = 1, 
h = 1, ft = 0.2 and n = 1000 for (a) 4 x 4 lattice and (b) 8 x 8 lattice, showing that 
Ph(M) = P eq (M), i.e., eigenfunction is indeed an equilibrium distribution. 



6. Summary 

In this paper we show and verify numerically that the repeated nonequilibrium 
measurements of work done to connect any two microstates of a system can be used to 
construct a matrix S whose principal eigenvector is the equilibrium distribution. The 
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matrix elements of S (Eq. [25]) for a Hamiltonian H(A,x) with (A,x) as a conjugate 



where the summation is over all paths that start from an equilibrium distribution of 
externally controlled parameter A = Ao with value of conjugate variable x as Xi and 
end in a state with A = A and x — Xf, with proper normalization. The work done 
W is defined in Eq. CD The values of the elements of S depend on the details of 
the process and, therefore, there can be many different S, but all will have the same 
invariant principal eigenvector. In this way the distribution of an interacting system 
can be obtained from a free, non-interacting one without any reference to equilibrium 
anywhere. In the process, we also provide a dynamics independent proof of the result 
that the equilibrium probability distribution can be obtained using the nonequilibrium 
path integral. Besides giving a new perspective of thermodynamics and statistical 
mechanics, our result has direct implications for new ways in numerical simulations 
and experiments. 
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